
12/26/13 End of day:

1) linear single vector is as fast as 3-array <array>

	--> how would usage work with single vector??

	--> or use separate x,y,z vectors??

2) <array> : fixed size!!!
	
	DONE: --> need new[] and del[] old-style array in SasAtom class

	--> try BOOST

	--> test both

3) need to finish dcd methods (sasio.py-like usage)

	--> read and write

	--> with accessor / step-wise methods

4) need to see how map : dictionary works for c++
	
	--> get mass

5) usage & inheritance (serial vector or array)
	
	--> calc_com

6) pdb read and write


	--> calc_rg

7) need to work on sassubset
 
	--> general slicing
	--> eval syntax




to find duplicate symbols

nm -j *.o | sort | uniq -d



### clean use of template ...

#include <iostream>
#include <vector>

using namespace std;

const int N = 40;

template <typename T>
int sum(T val)
{
	int s = 0;
	for(int i = 0; i < val.size(); ++i)
		s += val[i];
	return s;
}

int main()
{
	vector<int> vec;
	
	for(int i = 0; i < N; ++i)
		vec.push_back(i);
		
	cout<<" sum is "<< sum(vec) << '\n';
	return 0;
}
